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Abstract 

Charge collection measurements performed on heavily irradiated p-spray dofz pixel 
sensors with a grazing angle hadron beam provide a sensitive determination of 
the electric field within the detectors. The data are compared with a complete 
charge transport simulation of the sensor which includes signal trapping and charge 
induction effects. A linearly varying electric field based upon the standard picture 
of a constant type-inverted effective doping density is inconsistent with the data. 
A two-trap double junction model implemented in ISE TCAD software can be 
tuned to produce a doubly-peaked electric field which describes the data reasonably 
well at two different fluences. The modeled field differs somewhat from previous 
determinations based upon the transient current technique. The model can also 
account for the level of signal trapping observed in the data. 



1 Introduction 



In recent years, detectors consisting of one and two dimensional arrays of 
silicon diodes have come into widespread use as tracking detectors in particle 
and nuclear physics experiments. It is well understood that the intra-diode 
electric fields in these detectors vary linearly in depth reaching a maximum 
value at the p-n junction. The linear behavior is a consequence of a constant 
space charge density, iV eff , caused by thermodynamically ionized impurities 
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in the bulk material. It is well known that the detector characteristics are 
affected by radiation exposure, but it is generally assumed that the same 
picture is valid after irradiation. In fact, it is common to characterize the 
effects of irradiation in terms of a varying effective charge density. The use 
of N e g to characterize radiation damage has persisted despite a growing body 
of evidence [l]-[6] that the electric field does not vary linearly as a function 
of depth after heavy irradiation but instead exhibits maxima at both n+ and 
p+ implants. The work presented in this paper demonstrates conclusively 
that the standard picture does not provide a good description of irradiated 
silicon pixel detectors. We show that it is possible to adequately describe 
the charge collection characteristics of a heavily irradiated silicon detector 
in terms of a tuned double junction model which produces a doubly-peaked 
electric field across the detector. The allowed parameter space of the model 
can also accommodate the expected level of leakage current and the level of 
signal trapping actually observed in the detector. 

This paper is organized as follows: Section 2 describes the experimental tech- 
nique and data, Section 3 describes the carrier transport simulation used to 
interpret the data, Section 4 describes the technique used to model doubly- 
peaked electric fields and the limitations of previous models, Section 5 de- 
scribes the tuning of a successful model and describes its consequences, and 
Section 6 summarizes the results and develops several conclusions. 



2 Experimental Technique and Data 

This investigation is based upon beam test data that were accumulated as 
part of a program to develop a silicon pixel tracking detector [7] for the CMS 
experiment at the CERN Large Hadron Collider. The measurements were 
performed in the H2 line of the CERN SPS in June and September 2003 using 
150-225 GeV pions. The beam test apparatus is described in Ref. [9] and is 
shown in Fig. 1. A silicon beam telescope [8] consisted of four modules each 
containing two 300 /im thick single-sided silicon detectors with a strip pitch 
of 25 /im and readout pitch of 50 /jm. The two detectors in each module 
were oriented to measure horizontal and vertical impact coordinates. A pixel 
hybrid detector was mounted between the second and third telescope modules 
on a cooled rotating stage. A trigger signal was generated by a silicon PIN 
diode. The analog signals from all detectors were digitized in a VME-based 
readout system by two CAEN (V550) and one custom built FADCs. The 
entire assembly was located in an open-geometry 3T Helmholtz magnet that 
produced a magnetic field parallel to the beam. The temperature of the test 
sensors was controlled with a Peltier cooler that was capable of operating down 
to -30°C. The telescope information was used to reconstruct the trajectories 
of individual beam particles and to achieve a precise determination of the 
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particle hit position in the pixel detector. The resulting intrinsic resolution of 
the beam telescope was about 1 /im. 




side view 




Cooling system 



Fig. 1. The beam test apparatus consisting of four horizontal and four vertical planes 
of silicon strip detectors and a rotating stage for the pixel detector. 



2.1 Pixel Hybrids 



The prototype pixel sensors are so-called "n-in-n" devices: they are designed 
to collect charge from n+ structures implanted into n- bulk silicon. This de- 
sign is thought to provide greater longevity in high radiation fields and to 
allow "partially-depleted" operation after "type-inversion" of the substrate. It 
is more expensive than the "p-in-n" process commonly used in strip detectors 
because it requires double-sided processing and the implementation of inter- 
pixel isolation. Two isolation techniques were tested: p-spray, where a uniform 
medium dose of p-impurities covers the whole structured surface, and p-stop, 
where higher dose rings individually surround the n+ implants. Results on 
the Lorentz angle and charge collection efficiency measurements as well as 
a detailed description of both designs can be found elsewhere [9,10]. In this 
paper we discuss only measurements performed on p-spray sensors. All test 
devices were 22x32 arrays of 125x125 /xm 2 pixels having a sensitive area of 
2.75x4 mm 2 . The substrate was 285 jum thick, n-doped, diffusively-oxygenated 
silicon of orientation (111) and resistivity 2-5 kQ-cm. Individual sensors were 
diced from fully processed wafers after the deposition of under-bump metal- 
ization and indium bumps. A number of sensors were irradiated at the CERN 
PS with 24 GeV protons. The irradiation was performed without cooling or 
bias. The applied fluences 1 were 8.1xl0 14 n eq /cm 2 and 9.7xl0 14 n eq /cm 2 . In 
order to avoid reverse annealing, the sensors were stored at -20° C after irra- 
diation and kept at room temperature only for transport and bump bonding. 
All sensors were bump bonded to PSI30/AC30 readout chips [11] which allow 

1 All particle fluences are normalized to 1 MeV neutrons (n eq /cm 2 ). 
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analog readout of all 704 pixel cells without zero suppression. The PSI30 also 
has a linear response to input signals ranging from from zero to more than 
30,000 electrons. 



2.2 Data 



The main focus of the work presented in this paper involves a set of charge 
collection measurements that were performed using the "grazing angle tech- 
nique" [12]. As is shown in Fig. 2, the surface of the test sensor is oriented by 
a small angle (15°) with respect to the hadron beam. A large sample of data 
is collected with zero magnetic field. The charge measured by each pixel along 
the y direction samples a different depth z in the sensor. Precise entry point 
information from the beam telescope is used to produced finely binned charge 
collection profiles. For unirradiated sensors, the cluster length determines the 
depth over which charge is collected in the sensor. 



, z-axis 




y-ax is 

sensor backplane 

pixel implant 
bump bond 



Fig. 2. The grazing angle technique for determining charge collection profiles. The 
cluster length is proportional to the depth over which charge is collected. 



The profiles that were observed for an unirradiated sensor and for a sensor 
that was irradiated to a fluence of $ = 8.1 x 10 14 n eq /cm 2 are shown in 
Fig. 3. The unirradiated sensor was operated at a bias voltage of 150V which 
is well above its depletion voltage (approximately 70V). It produces a large 
and uniform collected charge distribution indicating that it is fully depleted (a 
large y coordinate indicates a large collection distance). The irradiated sensor 
was operated at bias voltages varying from 150V to 600V. It appears to be 
partly depleted at 150V, however, signal is collected across the entire thickness 
of the detector. Another puzzle is that the ratio of charges collected at 300V 
bias and 150V bias is Q (300V) /Q{ 150V) = 2.1 which is much larger than the 
maximum value of \/2 expected for a partly depleted junction. It is clear that 
the profiles for the irradiated sensor exhibit rather different behavior than one 
would expect for a heavily- doped, unirradiated sensor. 
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Fig. 3. Charge collection profiles for an irradiated ($=8x 10 n eq /cm 2 ) and an 
unirradiated sensor ($ = n cq /cm 2 ) operated at several bias voltages. 

3 Simulation and comparison with data 



It is well-known that signal trapping is a significant effect in heavily irradiated 
silicon detectors. In order to evaluate the effects of trapping, it is necessary to 
implement a detailed simulation of the sensor. Our simulation, Pixelav [13], 
incorporates the following elements: an accurate model of charge deposition 
by primary hadronic tracks (in particular to model delta rays) ; a realistic elec- 
tric field map resulting from the simultaneous solution of Poisson's Equation, 
carrier continuity equations, and various charge transport models; an estab- 
lished model of charge drift physics including mobilities, Hall Effect, and 3-d 
diffusion; a simulation of charge trapping and the signal induced from trapped 
charge; and a simulation of electronic noise, response, and threshold effects. 
A final step reformats the simulated data into test beam format so that it can 
processed by the test beam analysis chain. 

Several of the Pixelav details described in Ref. [13] have changed since they 
were published. The commercial semiconductor simulation code now used to 
generate a full three dimensional electric field map is the ISE TCAD package 
[14] . The charge transport simulation was modified to integrate only the fully- 
saturated drift velocity, 



dr fx[qE + fxr H E x B + qfi 2 r 2 H (E ■ B)B 
dt = ~ ^ 



l + /i 2 r 2 H \B\ 



where fi(E) is the mobility, q — ±1 is the sign of the charge carrier, E is 
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the electric field, B is the magnetic field, and r# is the Hall factor of the 
carrier. The use of the fully saturated drift velocity permits much larger inte- 
gration steps (which had previously been limited by stability considerations) 
and significantly increases the speed of the code. 

The simulation was checked by comparing simulated data with measured data 
from an unirradiated sensor. A plot of the charge measured in a single pixel 
as a function of the horizontal and vertical track impact position for normally 
incident tracks is shown in Fig. 4. The simulation is shown as the solid his- 
togram and the test beam data are shown as solid points. Note that the sensor 
simulation does not include the "bias dot" structure on the n+ implants which 
is used to provide a high resistance connection to ground and to provide a test 
pulse capability. There is reduced charge collection from this portion of the 
implant and the data shows reduced signal in both projections at the bias dot. 
Another check, shown in Table 1, is the comparison of the average Lorentz an- 
gle measured at several bias voltages [9] . In both cases, reasonable agreement 
is observed between measured and simulated data. 
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Fig. 4. A plot of the charge measured in a single pixel as a function of the hori- 
zontal and vertical track impact position for tracks that are normally incident on 
an unirradiated sensor. The simulation is shown as a solid histogram and the test 
beam data are shown as solid points. 

Table 1 

Measured and simulated values of average Lorentz angle 9l versus bias voltage for 
an unirradiated sensor. 



Bias Voltage 


Measured 9l [deg] 


Simulated 0l [deg] 


150V 


22.8 ±0.7 


24.7±0.9 


300V 


14.7 ±0.5 


17.4±0.9 


450V 


11.2 ±0.5 


12.0±0.9 



6 



The charge collection profiles for a sensor irradiated to a ffuence of $ = 
8.1 x 10 14 n eq /cm 2 and operated at bias voltages of 150V and 300V are 
presented in Fig 5. The measured profiles are shown as solid dots and the 
simulated profiles are shown as histograms. The simulated profiles were gen- 
erated with electric field maps that correspond to two different effective den- 
sities of acceptor impurities. The red histograms are the simulated profile for 
N c g = 3.5 x 10 12 cm -3 . Note that the 300V simulation agrees well with the 
measured profile but thel50V simulation is far too broad. The blue histograms 
show the result of increasing N e s to 24 x 10 12 cm -3 . At this effective doping 
density, the width of the simulated peak in the 150V distribution is close to 
correct but it does not reproduce the "tail" seen in the data at large y. The 
300V simulated distribution is far too narrow and the measured signals are 
too small (note that the profiles are absolutely normalized). It is clear that a 
simulation based upon the standard picture of a constant density of ionized 
acceptor impurities cannot reproduce the measured profiles. 




Position ((.im) 



Position ((.im) 



Fig. 5. The measured and simulated charge collection profiles for a sensor irradi- 
ated to a ffuence of $ = 8.1 X 10 14 n eq /cm 2 . The profiles measured at bias voltages 
of 150V and 300V are shown as solid dots. The red histogram are the simulated 
profiles for a constant effect doping N e fi = 3.5 x 10 12 cm~ 3 of acceptor impuri- 
ties. The blue histograms are the simulated profiles for a constant effect doping 



N, 



eff 



24 x 10 12 cm 



Note that the simulation of the irradiated sensor includes the effects of trap- 
ping. The trapping rates of electron and holes have been shown to scale linearly 
with ffuence [15,16]. Unfortunately, the fluence measurements have a fractional 
uncertainty of ±10% which feeds directly into an uncertainty on the trapping 
rates. An additional uncertainty arises because annealing can modify the trap- 
ping rates by 30% [15] leading to a fairly large overall uncertainty. The charge 
collection profiles are not very sensitive to the hole trapping rates but are 
sensitive to the electron trapping rates. The profiles shown in Fig. 5 are simu- 
lated by reducing the trapping rates by 20% with respect to the nominal rates 
for reasons that will be explained in Section 5. The choice does not affect the 
conclusion that the measured profiles are inconsistent with a constant effective 
doping model. 
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4 Double Junction Models 



The large number of measurements suggesting that large electric fields exist 
on both sides of an irradiated silicon diode has given rise to several attempts 
to model the effect [3,6,17]. The most recent of these by Eremin, Verbitskaya, 
and Li (EVL) [17] is based upon a modfication of the Shockley- Read-Hall 
(SRH) statistics. The EVL model produces an effective space charge density 
Pes from the trapping of leakage current by one acceptor trap and one donor 
trap. The effective charge density is related to the occupancies and densities 
of traps as follows, 

Pcff = e [N D f D - N A f A ] + Pdopants (2) 



where: N D and N A are the densities of donor and acceptor trapping states, 
respectively; fr> and f A are the occupied fractions of the donor and acceptor 
states, respectively; and Pdopants is the charge density due to ionized dopants. 
Charge flows to and from the trapping states due to generation and recom- 
bination. The occupied fractions are given by the following standard SRH 
expressions, 



v h a%p + v e a D 'n i e ED/kT 
v e a E (n + me E °/ kT ) + v h a%{p + me- E °l kT ) 

Ve^n + v h a^nie- EA/kT 

v e a^(n + me EA l kT ) + v h afi(p + me- EA / kT ) 



where: v e and Vh are the thermal speeds of electrons and holes, respectively; 
a® and a E are the electron and hole capture cross sections for the donor trap; 

and are the electron and hole capture cross sections for the acceptor 
trap; n and p are the densities of free electrons and holes, respectively; rii is the 
intrinsic density of carriers; E D and E A are the activation energies (relative 
to the mid-gap energy) of the donor and acceptor states, respectively. The 
generation-recombination current caused by the SRH statistics for single donor 
and acceptor states is given by the following expression, 



VhVeOh°®N D {rip - nf) 

v e a^(n + me ED l kT ) + v h a^(p + me- ED l kT ) 

v h Ve(Th a t N A(np - nf) 

v e a^(n + nie EA / kT ) + v h a£(p + riie- EA l kT ) ' 



Within the EVL model, the four trapping cross sections are fixed to 10 15 cm 2 . 
The leakage current is generated from an additional SRH trapping state that is 
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introduced for this purpose but is assumed not to trap charge. The donor and 
acceptor states are assumed not to generate leakage current which, given the 
small size of the cross sections, is a self-consistent assumption. The densities 
of the donor and acceptor states (N D and N A ) are adjusted to "fit" TCT data. 
The parameters of the model are given in Table 2. The trap densities are scaled 
to fluence and are given in terms of introduction rates g 1Ilt = NA/o/^eq- 

Table 2 

Parameters of the EVL Model [17]. 



Trap 


E (eV) 


ffint (cm" 1 ) 


a e (cm 2 ) 


a h (cm 2 ) 


Donor 


E v + 0.48 


6 


1 x 10~ 15 


1 x 10~ 15 


Acceptor 


E c - 0.525 


3.7 


1 x 10~ 15 


1 x 10~ 15 



An illustrative sketch of the EVL model has been reproduced from Ref. [17] 
and is shown in Fig. 6. Figure 6a shows a uniform current density flowing across 
a reverse-biased junction. Since holes are produced uniformly across the junc- 
tion and flow to the p+ backplane, the hole current density increases linearly 
with increasing z from the n+ implant to the p+ implant. The electrons flow 
to the n+ implant and the electron current density increases with decreasing 
z. The actual carrier densities depend upon the details of the fields and mo- 
bilities but vary monotonically across the junctions as shown in Fig. 6b. The 
trapping of the mobile carriers produces a net positive space charge density 
near the p+ backplane and a net negative space charge density near the n+ 
implant as shown in Fig. 6c. Since positive space charge corresponds to n-type 
doping and negative space charge corresponds to p-type doping, there are p-n 
junctions at both sides of the detector. The electric field in the sensor follows 
from a simultaneous solution of Poisson's equation and the carrier continuity 
equations. The resulting z-component of the electric field is shown in Fig. 6d. 
It varies with an approximately quadratic dependence upon z having a min- 
imum at the zero of the space charge density and maxima at both implants. 

In order to test whether the electric field predicted by the EVL model would 
improve the agreement between the simulated and measured charge collection 
profiles shown in Section 3, it was necessary to implement the EVL model 
in ISE TCAD. TCAD contains a complete implementation of SRH statistics. 
However, the EVL modifications of SRH are not incorporated. In particular, 
any state added to generate leakage current would also trap charge. It is 
possible to use the TCAD "Physical Model" interface to replace the existing 
SRH implementation with a modified one, however, a less invasive approach 
was adopted. 

Our approach is based upon three observations: 
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Fig. 6. An illustrative sketch of the EVL model [17]. 



(1) The trapping cross sections are poorly known. The cross sections for 
states observed in various types of defect spectroscopy vary over several 
orders of magnitude. 

(2) The occupancies fo/A of the trapping states are independent of the scale 
of the cross sections. If the capture cross sections a e /h in Eq. 3 are rescaled 
by a factor r, then fn/A are unchanged. The occupancies depend only 
upon the ratio of electron and hole cross sections, cr^/ce. 

(3) The current U generated by the donor and acceptor impurities is linear 
in the cross section rescaling factor r. 

These observations imply that it is possible to implement the EVL model in 
TCAD simply by setting a® = a® = = and by varying the size of the 
common cross section until the generation current is equal to the observed 
or expected leakage current. The trap occupancies are not affected in zeroth 
order by the rescaling, but the leakage current and the free carrier densities are 
affected by r. The carrier densities have a first-order effect on the occupancies 
so that varying r does alter p e g. This approach uses the same trapping states 
to produce space charge and leakage current (it is not necessary to introduce 
current-generating states). 

Using the activation energies and introduction rates for the donor and acceptor 
states given in Table 2, the only free parameter in the TCAD implementation 
of the EVL model is the size of the common cross section or equivalently, 
the leakage current. The leakage current in the test sensors was observed to 
increase substantially after bump-bonding to the readout chips. Although the 
cause of the increased current is not clear, it is possible that it is due to 
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increased surface/edge leakage or that it is caused by stressing the detector. 
The measured leakage current, expressed in terms of the damage parameter a 
[18], is listed in Table 3. It is well established [18] that in large detectors, the 
magnitude of leakage current is expected to be o — 4 x 1CT 17 A/cm for a fully 
biased detector. It is clear that the observed leakage current is approximately 
five times larger than the expected value. 

Table 3 

Observed leakage current in terms of a = li ca k(20 o C)/(V3>) where: /i C ak(20°C) is 
the leakage current expected at temperature 20°C, V is the volume of the sensor, 
and $ is the neutron-equivalent fluence. 



Bias Voltage 


a 


(A/cm) 


-150V 


15 


x 10" 


-17 


-300V 


19 


x 10" 


-17 


-450V 


25 


x 10" 


-17 



The EVL model was implemented in the sensor simulation using several pa- 
rameter choices. The resulting charge collection profiles for 150V and 300V 
bias voltages are shown in Fig. 7. The measured profiles are again shown 
as solid dots. The blue histogram shows the EVL model with the leakage 
current adjusted to 20% of the measured value which is comparable to the 
expected leakage current (cto). It clearly substantially underestimates the col- 
lected charge at both voltages. The effect of increasing the leakage current to 
60% of the measured value is shown as the solid red histogram. This improves 
the agreement at 300V but produces too much signal in the high z tail of 
the 150V distribution. Finally, in an attempt to increase the collected charge 
at 300V, the introduction rates were scaled down by a factor of eight. The 
common cross section was increased by the same factor to hold the leakage 
current constant. The result of this is shown as the black histogram in Fig. 7. 
Note that this does increase the collected charge signal, however, it is still too 
small at 300V and is much too large in the tail region of the 150V distribution. 
We conclude that the EVL model does not provide a quantitative description 
of the measured charge collection profiles. 



5 An Improved Two- Trap Model 

Although it led to a poor quantitative description of the measured charge col- 
lection profiles, the EVL model does have the qualitative features needed to 
describe the data. At low bias voltages, the combination of the quadratic min- 
imum in the electric field and signal trapping can act like a "gate" suppressing 
the collection of charge from the p+ side of the detector. The measured profile 
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Fig. 7. The measured charge collection profiles (solid dots) at 150V and 300V are 
compared with simulations based upon the EVL model. The blue histogram shows 
the EVL model with the leakage current adjusted to 20% of the measured value 
which is comparable to the expected leakage current («o)- The red histogram shows 
the effect of increasing the leakage current to 60% of the measured value. The black 
histogram shows the result of decreasing the introduction rates by a factor of eight 
and increasing the cross sections by the same factor. 



would then appear to be a narrow peak on the n+ side of the detector. As 
the bias is increased, the magnitude of the field at the minimum would also 
increase and effectively "lift the gate" which would allow much more charge 
collection from the p+ side of the detector. 

In order to investigate whether a two-trap double junction EVL-like model 
can describe the measured charge collection profiles, a tuning procedure was 
adopted. Relaxing the EVL requirement that all trapping cross sections are 
equal, the model has six free parameters (No, Na, <J®, <Jh ', u^, (?h) that can 
be adjusted. The activation energies are kept fixed to the EVL values. Addi- 
tionally, as was discussed in Section 3, the electron and hole signal trapping 
rates, T e and T h , are uncertain at the 30% level due to the fluence uncertainty 
and possible annealing of the sensors. They are treated as constrained param- 
eters. In an ideal case, the simulation would be used to perform a x 2 fit to the 
measured profiles and optimal parameter estimates and confidence intervals 
would be extracted. Unfortunately, the simulation chain is quite slow. It re- 
quires approximately nine hours to produce a TCAD solution for the electric 
field map and then several simultaneous Pixelav jobs (one per bias voltage) re- 
quire approximately 16 hours each to produce adequate statistics. This implies 
that a traditional x 2 fit with a full error matrix analysis is not feasible. An 
extremely tedious and less mathematically rigorous procedure was adopted in- 
stead. The parameters of the double junction model were systematically varied 
and the agreement between measured and simulated charge collection profiles 
was judged subjectively. 

In the course of the tuning procedure, it became clear that the EVL model does 
not produce a sufficiently large electric field on the p+ side of the detector. 
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The solution to this problem is to increase the density of donors (hole traps) 
as compared with the density of acceptors (electron traps). When this is done, 
the z position of the minimum in the effective charge density shifts toward the 
n+ implant as sketched in Fig. 8a. Unfortunately, this causes the "peak" in 
the 150V simulated charge profile to become too narrow. The position of the 
charge density minimum can be restored to a position nearer the midplane of 
the detector by decreasing the ratios of the hole and electron cross sections 
from 1.0 to 0.3 (crf/cxf = 0.3 and cr^/a^ = 0.3) as shown schematically 
in Fig. 8b. Note that the adjustment of the cross sections for both types of 
trap minimizes the field in the quadratic minimum while allowing for large 
fields at the implants. For simplicity it is assumed that the electron cross 
sections are equal (af = = a e ) and that the hole cross sections are equal 



Fig. 8. The effect of increasing Nd/Na when (a) the electron and hole cross sections 
are equal, and when (b) ah/a e = 0.3. 

The current best "fit" to the measured charge collection profiles is called model 
dj35 and reduces the ratio Na/Nd from the EVL value of 0.62 to 0.4. The 
^-component of the simulated electric field is plotted as a function of z in 
Fig. 9 for bias voltages of 150V and 300V. The field profiles have minima 
near the midplane of the detector. Note that the minimum field at 150V bias 
appears to be very small but is still approximately 400 V/cm. The dj35 field 
profiles are compared with those that result from the rescaled EVL model 
which are shown as blue solid and dashed curves (labeled as dj 16) . Finally, 
the electric field profiles resulting from a constant p-type doping of density 
iVcfj = 3.5 x 10 12 cm -3 are shown as green curves for reference. 

The measured charge collection profiles at bias voltages of 150V, 300V, and 
450V are compared with the dj35 simulation in Fig. 10 for fluences of: (a) 
8.1 x 10 14 n eq /cm 2 and (b) 9.7 x 10 14 n eq /cm 2 . The electron trapping rate for 
the lower fluence is set to 80% of the nominal value. The electron trapping 
rate for the higher fluence sample which was held at 30° C for a longer time 
is set to 70% of the nominal value. The trap densities for the higher fluence 
simulation were scaled linearly from the densities used for the lower fluence 
simulation. 
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Fig. 9. The z-component of the simulated electric field resulting from model dj35 
is plotted as a function of z. The field profiles for 150V and 300V are shown as 
as red solid and dashed curves, respectively. The field profiles that result from the 
rescaled EVL model are shown as blue solid and dashed curves and labeled as 
djl6. The electric field profiles resulting from a constant p-type doping of density 
N e s = 3.5 x 10 12 cm~ 3 are shown as green curves. 

Although the 8.1 x 10 14 n eq /cm 2 fluence simulation falls a below the measured 
profile at small y (near the n+ implant), it provides a reasonable description of 
the measurements. Several features of the measured distributions are described 
well by the simulation. Note that both data and simulation show a distinctly 
negative signal near y — /im. This can be understood as a consequence 
of hole trapping. Electrons deposited near the n+ implant are collected with 
high efficiency whereas holes deposited near the implant must transit the entire 
detector thickness to reach the p+ backplane. If the holes are collected, they 
produce no net signal on the n+ side of the detector. However if the holes 
trap, then a negative signal is induced and is most visible in the y < jum 
region. Another feature is the "wiggle" in the 150V profiles for both low and 
high fluences. The relative signal minimum near y = 700 /im corresponds to 
the E z minimum where both electrons and holes travel only short distances 
before trapping. This small separation induces a only a small signal on the n+ 
side of the detector. At larger values of y, E z increases causing the electrons 
drift back into the minimum where they are likely to trap. However, the holes 
drift into the higher field region near the p+ implant and are more likely to 
be collected. The net induced signal on the n+ side of the detector therefore 
increases and creates the local maximum seen near y = 900 /mi. 

At the 9.7 x 10 14 n eq /cm 2 fluence, the simulation provides a poorer descrip- 
tion of the corresponding measurements than it does at the lower fluence. 
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Fig. 10. The measured charge collection profiles at bias voltages of 150V, 300V, 
and 450V are shown as solid dots for fluences of: (a) 8.1 x 10 n eq /cm 2 and (b) 
9.7 x 10 14 n eq /cm 2 . The dj35 simulation is shown as the red histogram in each plot. 

It is likely that there was more annealing of the higher fluence sample and 
simple linear scaling of trap densities with fluence is probably not optimal. 
The accommodation of the additional annealing may require an adjustment 
of the Na/Nd ratio or other degrees of freedom. A second observation is that 
the measured and simulated signals are absolutely normalized and the signal 
observed at 450V from the high fluence sample is larger than that observed 
for the lower fluence sample. This is surprising and may suggest that sample- 
to-sample normalization issues exist or it may confirm the idea that annealing 
has altered the behavior of the sample. 



The dj35 model fixes the ratios Na/N d and <7h/<y e leaving the parameters 
N D and a e to vary. Over a restricted range, increases in N D can be offset 
by decreases in a e keeping the electric field profile more or less unchanged. 
Scaling the electron cross section as a e oc N^ 2 ' 5 produces very similar charge 
collection profiles. The allowed region in N£,-a e space is shown in Fig. 11a as 
the solid line in the logarithmic space. If the donor density becomes too small 
(Nd < 20 x 10 14 cm -3 ), the 150V simulation produces too much signal at large 
z. If the donor density becomes too large (N D > 50 x 10 14 cm -3 ), the 300V 
simulation produces insufficient signal at large z. Since the simulated leakage 
current varies as Ji ea k oc u e N£,, different points on the allowed solid contour 
correspond to different leakage current. Contours of constant leakage current 
are shown as dashed curves and are labeled in terms of the corresponding 
damage parameter a where a is the expected leakage current. It is clear that 
the simulation can accommodate the expected leakage current which is smaller 



15 



than the measured current by a factor of 4-5. 




Fig. 11. The allowed region in Nu-a e space for model dj35 is shown as the solid 
red line in (a) and (b). Contours of constant leakage current are shown as dashed 
curves in (a) and are labeled in terms of the corresponding damage parameter a 
where ao is the expected leakage current. Contours of constant electron trapping 
rate are shown as dashed curves in (b) and are labeled in terms of the un-annealed 
trapping rate Tq for the nominal fluence. 



The electron and hole traps in the model should also contribute to the trapping 
of signal carriers. The contributions of these states to the effective trapping 
rates of electrons and holes are given by the following expression, 



■v e 
v h 



a°N D (l 



1 - f A ) + <r?N D f D } ~ v e a*N A 
f D )+a£N A f A 



(5) 



v h a%N D 



where it has been assumed that the trap occupancies are small. Because 
Na/Nd is assumed to be constant, contours of constant electron trapping 
rate are parallel to contours of constant leakage current in Nr,-a e space. The 
best "fit" of the simulation to the measured profiles reduced T e to 80% of the 
un-annealed trapping rate To for the nominal fluence [15] . These contours are 
compared with the allowed contour in Fig. lib. It is clear that the simulation 
can accommodate the measured trapping rate in the same region of parameter 
space that maximizes the leakage current. 

Figure lib also suggests a solution to the puzzle that the trapping rates have 
been shown to be unaffected by the presence of oxygen in the detector bulk 
[15] whereas it is well-established that the space charge effects are quite sen- 
sitive to the presence of oxygen in the material [19]. It is clear that additional 
small-cross-section trapping states can play a large role in the effective charge 
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density but a small one in the effective trapping rates. If the formation of 
the additional small-cross-section states were suppressed by oxygen, then p e ^ 
could be sensitive to oxygenation whereas T e / h would be insensitive to oxy- 
genation. This is another consequence of the observation that the occupancies 
}'d/a of the trapping states are independent of the scale of the cross sections 
in the steady state (see Section 4). The trapping of signal carriers is not a 
steady-state phenomenon and is sensitive to the scale of the trapping cross 
sections. 



6 Conclusions 

The main result of the work presented in this paper is that a doubly-peaked 
electric field is necessary to describe to the charge collection profiles measured 
in heavily irradiated pixel sensors . A simulation utilizing a linearly varying 
electric field based upon the standard picture of a constant type-inverted ef- 
fective charge density is inconsistent with the measurements. 

A two-trap EVL-like model can be tuned to provide a reasonable description 
of the measurements. It can also account for the expected level of the leakage 
current (although not the observed level of leakage current) and the observed 
electron signal trapping rate. It is important to state that any two-trap model 
is, at best, an "effective theory". It is well-known that there are a great many 
trapping states in heavily irradiated silicon. All of them can and do trap 
charge. There may also be thermodynamically ionized defects that contribute 
to the effective space charge density. It is clear that a two-trap model can 
describe the gross features of the physical processes but it may not be able 
to describe all details. This also implies that the parameters of the two-trap 
model presented in this paper are unlikely to have physical reality. 

The charge-sharing behavior and resolution functions of many detectors are 
sensitive to the details of the internal electric field. A known response function 
is a key element of any reconstruction or analysis procedure. A working effec- 
tive model will permit the detailed response of these detectors to be tracked 
as they are irradiated in the next generation of accelerators. 

Finally, we note that despite a growing body of contrary evidence, the overly 
simple picture of uniform type inversion in irradiated silicon detectors per- 
sists in the minds of many researchers and even in the terminology used by 
experts to describe their devices. After irradiation, quantities like £dep (de- 
pletion depth) and N e g may correctly suggest reduced detector performance 
but given the reality of doubly-peaked electric fields and signal trapping, they 
have no physical significance. 
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